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Abstract 

The renormalization, at the next-to- leading order in a s , of the AB = 2 operators at the 
^ ■ lowest order in the heavy quark expansion, namely in the static theory, is computed taking 
into account previously missed contributions. These operators are relevant for the calculation 
of the B°-B° mixing on the lattice. 
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1 Introduction 



The lack of a precise knowledge of the matrix elements of the operator 

(B\h»{l - l5 )d fc 7/i (l - % )d\B) = p 2 B B B M B , (1) 

generated by the box diagram with the exchange of virtual top quarks, is the main source of 
theoretical error in the extraction of the CKM matrix element Vtd from the B°-B° mixing pa- 
rameter Xd- This uncertainty on f B B B propagates to other theoretical estimates. In particular 
the extraction of the CP- violating phase from the combined analysis of Xd and the K°-K° CP- 
violating parameter e is affected by a twofold ambiguity that would be eliminated by a precise 
determination of f B B B . In turn this would strongly reduce the uncertainty in the prediction 
of the asymmetry in the decay B —>■ J ftp K$ [0. 

The matrix element in eq. ([I]) can be evaluated on the lattice [[|. However, since m?, is 
larger than the current values of the lattice cutoff, the h quark cannot be put on the lattice as a 
dynamical field. Therefore an effective theory based on the expansion in the heavy quark mass 
is needed. Such a theory has been built 0, and its discretized version can be used in lattice 
simulations. In particular the expansion at the lowest order in is used to build the static 
effective theory both in the continuum and on the lattice. This theory has the usual form 

7~tstatic — ^ CiQi (2) 

i 

in terms of Wilson coefficients and local four-fermion operators. In order to use lattice results, 
the effective theory in the continuum must be matched both to its lattice counterpart and to 
the "full" theory, namely a theory with a dynamical heavy field. In ref. |4]] this matching has 
been computed at 0(a s ), using the AB = 2 effective Hamiltonian || as the "full" theory. 

The two matching procedures take place at different scales. In fact the matching to the "full" 
theory is done at a scale of the order of the ultraviolet cutoff, namely m^, while the matching 
to the static theory on the lattice is done at typical current values of the lattice cutoff, 1/a ~ 2 
GeV. A complete determination of the static theory in the continuum requires the calculation 
of the running of the Wilson coefficients between these two scales. This is usually done by 
using the renormalization group equations (RGEs). In this way, at the leading order (LO), 
one resums in the Wilson coefficients terms such as a™log n , assumed to be of 0(1). To be 
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consistent, an 0(a s ) matching calls for a next-to-leading (NLO) determination of the Wilson 
coefficients, which resums terms of the type log™. Of course, in the case at hand, one 
may argue that the relevant logs, namely log(a 2 m 2 ) ~ 1.6, are not large, so that the running 
at the leading order can be considered as a pure 0(a s ) effect. Anyway the calculation of the 
anomalous dimension up to the NLO is required to have a regularization-scheme-independent 
expression of the Wilson coefficients. In the past an effort has been made to calculate this 
NLO anomalous dimension in the static theory ||. However, as pointed out in ref. f7|, some 
contributions coming from the operator mixing have been overlooked in ref. . 

In the next section we calculate these new contributions and present a complete determi- 
nation of the NLO Wilson coefficients of the static theory. 

2 NLO Wilson coefficients in the static theory 

In this section we discuss the NLO gluon renormalization of the AB = 2 operators at the lowest 
order in the heavy quark expansion, i.e. in the static limit m& — > oo. We want to calculate the 
NLO expression of the Wilson coefficients of the relevant operators. To this end a few steps 
are required. First of all, the basis of local operators in the effective theory must be identified. 
Then the effective theory has to be matched against a "full" theory at a scale of the order 
of the ultraviolet cutoff, fixing in this way the initial conditions of the renormalization group 
equations. Finally the anomalous dimension matrix in the effective theory must be calculated 
at the desired order in a s and the RGE solved to give the Wilson coefficients as functions of the 
renormalization scale /i 2 . The scale- independent physical amplitude is given by the product 
of the Wilson coefficients C(fi 2 ) and the matrix elements of the corresponding renormalized 
operators (Q(fi 2 )), the latter usually requiring some non-perturbative method to be evaluated. 

The AB = 2 operator basis in the static limit is given by two dimension-six local operators 



Qi = 2/i (+) 7m (1 - l?,)d /i<-y (i - l5 )d , Q 2 = 2h™(l - 75 )rf U~\l - j s )d , (3) 



lation explicitly shows that indeed the basis is closed under renormalization at the NLO. 




where the field creates a heavy quark and •* annihilates a heavy antiquark. Our calcu 
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The next step is the NLO matching. The heavy quark theory has to be matched against 
the AB = 2 effective Hamiltonian by comparing the matrix elements at the scale mb of the 
relevant operators in the "full" and the effective theories, up to and including 0(a s ) terms. 
Here the effective Hamiltonian plays the role of the "full" theory, even if it is also an effective 
theory with the top and the heavy bosons integrated out, which in turn is matched against the 
Standard Model at the weak scale. The AB = 2 effective Hamiltonian has been calculated in 
ref. [|[ and is completely known at the NLO. There exists only one AB = 2 operator in this 
theory, namely 

QLL = h„(i-i 5 )dh>*(i-i 5 )d. (4) 

The calculation of the matching of this operator onto the operators in eq. @ requires the 
expansion of the matrix element (bd\Q LL\bd) in the heavy quark mass. This calculation has 
been done in refs. M, M and gives the initial condition at the scale mb 
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-14 NDR 

Bl ~\ -11 DRED ' ^--8. (5) 

We have reported the values of B\ both in the naive dimensional regularization scheme (NDR) 
and in dimensional reduction (DRED). This scheme dependence cancels out against the cor- 
responding dependence contained in the Wilson coefficient C^irn^) of the AB = 2 effective 
Hamiltonian^. The other initial condition, C^r??!), is the same both in DRED and NDR. We 
notice that, starting at 0(a s ), it does not pick up the scheme-dependent part inside Cxx,(m|) 
at the NLO. 

The evolution of the Wilson coefficients between the matching scale m% and the renormal- 
ization scale fi 2 is determined by the renormalization group equation 

A^<V) = ^<V) . (6) 
The anomalous dimension matrix 7 is defined as 

7 = 2A* 2 ^Z , (7) 

: As noticed in ref. 0, this implies that the anomalous dimension matrix element 7^, see eq. (pi]), in the 
effective theory is the same in DRED and NDR. 
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Figure 1: The Feynman diagrams that contribute to Thick (thin) lines represent the heavy 
(light) quarks. The blobs are the operator insertion point. 



where Z is the matrix of the renormalization constants connecting the bare and the renormalized 
operators 

6 R = Z-'Ob . (8) 



The formal solution of eq. (0) is 



C(fi 2 ) = T as e 



Ja s (r. 



(9) 



where T as is the ordered product with increasing powers of the coupling constant from left to 
right and (3{a s ) is the QCD beta function 



(3{a s ) = ^ 



dfi 2 



(10) 



In order to calculate the Wilson coefficients at the NLO, the first two terms of the pertur- 
bative expansion of (3{a s ) and j(a s ) are needed: 
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Table 1: Po/e coefficients of the operator insertions into the diagrams of fig. 0, calculated in 
the linear gauge. An overall factor a s / 'An is understood. The second column contains the 
diagram multiplicity factors, which have already been applied to the shown coefficients. To 
account for the renormalization of the external legs, self-energy diagrams count as 1/2. 



7(a.) = ?7 (0) +f^ 2 7 (1) + ... • (11) 
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The beta function coefficients are well known: 

HJV-2n/ 34 „ 10 

A> = 3 — f - , Pi = -j N - Y Nnf ~ 2Cpnf ' (12) 

where N is the number of colours, nj is the number of active flavours and Cp = {N 2 — 1)/2N. 
The anomalous dimension matrices 7^ and are calculated by computing the one- and 
two-loop renormalization of the operators in eq. (||). 

The one-loop renormalization matrix is given by the infinite parts of the operator insertion 
into the diagrams in fig. Ill, computed in the static theory. We have used the dimensional 
regularizationQ to calculate the divergent parts of these diagrams that appear as poles in e = (4— 
D)/2. The coefficients of the poles are collected in table [I]. The calculation is straightforward, 
the only peculiarity being that few tensor structures can appear in the effective theory because 
of the equation satisfied by the static fields, v^^hy 1 ' = ±h^\ In particular tensors can be 
reduced in the following way: 

fc< + V(l- 7fi )d/^<7' B '(l-75)d = 4p + V(l- 75 )d^7/,(l- 7 5)d 

+^(+)(l- 75 )d^(-)(l- 75 )d] . (13) 

The anomalous dimension at the LO is minus twice the pole coefficients in the matrix Z of 



75 and subtraction prescriptions are immaterial for the anomalous dimension at the LO. 
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the operator renormalization constants, so that, from table [I], we readily find 
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(o) 



-6C7 



1 + 



N 





-6CV + 4 + 



N 



(14) 



As expected, the anomalous dimension matrix is independent of the gauge. 7^ agrees with 
the previous calculations M, while the other matrix elements were not considered previously. 
The form of this matrix has some important consequences. On the one hand, Q\ insertion has 
vanishing component on Q2 at the leading order. Since the initial condition of Q2 is already of 
order a s , the one-loop contribution of Q2 becomes a NLO effect. For the very same reason, we 
do not need to calculate the two-loop renormalization, i.e. the second row of j^, which would 
generate a next-to-next-to-leading order term. On the other hand, Q2 has a non-zero leading 
component on Q\, and thus contributes to the Wilson coefficient C\ at the NLO. 

We still have to calculate the first row of This is a hard task, involving the evaluation 
of the pole parts of several two-loop diagrams. However the renormalization of Q\ onto itself 
has already been calculated in ref. ||, while the insertion of Q\ has vanishing component onto 
Q2- This last statement holds to all orders in perturbation theory, provided that one chooses a 
renormalization scheme that preserves the Fierz symmetry. In fact 



Q (+) = Qi 



Q H = Q2 + -Qi 



(15) 



are the eigenvectors of the Fierz transformation with eigenvalues ±1 0, therefore they cannot mix 
under renormalization. This enforces the following relations among the anomalous dimension 
matrix elements 

712 = , 721 = -(722 - 7n) , (16) 
which indeed are satisfied by eq. (|14|). 

The NLO anomalous dimension matrix is then given by 
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AT-l 
' 127V 



127N 2 + U3N + 63 - f + 8vr 2 (n 2 - 2N + ±) - n f (28N + 44) 



7« 
7n 



where, as already noticed, the entry marked with X is not needed at the NLO. 
3 This is a consequence of the properties of the static fields, v^^h^ = ±.h^ . 





X 



(17) 
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We are now ready to write down the solution of the RGE, eq. (^), at the NLO. Using 
eqs. ([ID, (HD and (0), we obtain 



c * 2) = fig) 




where 




(19) 



The new contribution to the NLO running of Ci is the term proportional to 7^ , while C2 is 
renormalized multiplicatively. 

Numerically the new term contributes to the running between mf and a typical lattice scale 



enhancement coming from 7^ . Moreover the operator Q2 also contributes at the NLO and 
should be included in lattice calculations of j\B^ in the static limit. 
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